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ON THE NATURE OF THE FAST MOVING STAR S2 IN THE GALACTIC CENTER 
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ABSTRACT 

We analyze the properties of the star S2 orbiting the supermassive black hole at the center of 
the Galaxy. A high quality SINFONI H and K band spectrum obtained from coadding 23.5 hours 
of observation between 2004 and 2007 reveals that S2 is an early B dwarf (B0~2.5V). Using model 
atmospheres, we constrain its stellar and wind properties. We show that S2 is a genuine massive 
star, and not the core of a stripped giant star as sometimes speculated to resolve the problem of star 
formation so close to the supermassive black hole. We give an upper limit on its mass loss rate, and 
show that it is He enriched, possibly because of the presence of a magnetic field. 
Subject headings: Stars: early-type — stars: fundamental parameters — Galaxy: center 



1. INTRODUCTION 

The central parsec of our Galax y hosts a larg e pop- 
ulation of young in assive stars (jKrabbe et al.1 Il995t 
iPaumard et al.l[2006f ). Their presence is puzzling since 
according to standard theories, star formation should not 
happen so close to the supermassive black hole SgrA*: 
the tidal forces are so large that any molecular cloud 
should be disrupted before being able to collapse (Morria 
Il993f l. For most of the young stars, the solution seems 
to be star formation in dense accretion disks as wit- 
nessed by the t wo counter-rotating flat stellar structures 
around SgrA* (iLev in fc Bclobor odovll200l iGenzel et all 
120031 : [Paumard et al.. .2006. ) . However, this attractive 
solution does not explain the so-called "S stars" , the 
group of objects located within one arcsecond (=0.038 
c) of the black hole. T heir orbits are randomly oriented 
Eisenhauer et al.l l2005l and thus do not fit the accre- 
tion disk scenario . Spectroscopically iden tified as B stars 
(|Ghez et al.ll2003HEisenhauer et aLll2005l ). their true na- 
ture as massive stars has been challenged: they have been 
proposed to be the hot, luminous cores of evolved red gi- 
ant and/or AGB stars the envelope of whic h has been 
stripped by tidal interact ions with SgrA* (jAlexanderJ 
120051: iDavies fc Kindl2005[) . In that case, the problem 
of their formation process vanishes since they might very 
well have formed away from the hostile environment of 
the black hole before having migrated towards the Galac- 
tic Center. Such a process is too slow to explain the pres- 
ence of short lived massive stars close to Sgr A*. Hence, 
these S stars constitute a "paradox of youth" (jGhez et al.l 
[200l . 

In this letter, we analyze the properties of the brightest 
S-star - S2 - and show that it is a genuine massive star. 



2. OBSERVATIONAL DATA AND MODELS 

To constrain the physical parameters of S2, a high 
quality spectrum is needed. Since the observation of 
the S stars is time consuming, the S/N ratio obtained 
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in one night is usually limited to 10 at maximum. To 
increase this ratio, we have co-add ed all the S2 sp ec- 
tra obtained with SINFONI (Eisenh auer et al.l 120031 ) in 
adaptive optics mode (0.0125" x 0.025") on the VLT since 
2004. This corresponds to a total integration time of 23.5 
hours in K band. Each spectrum was carefully extracted 
by selection of source and background pixels from which 
nebular contamination was removed. Note however that 
this is not such a critical problem for S2 which has a 
large radial velocity so that the position of the stellar 
absorption core is blueshifted beyond the nebular emis- 
sion. We have determined the combined spectrum by a 
cross-correlation technique. First, the position of the Br7 
line was estimated by a simple Gaussian fit to the line. 
From the such obtained radial velocities a first combined 
spectrum was calulated. This spectrum was then used 
to cross-correlate all individual spectra against it. This 
yielded better estimates for the radial velocities, which 
in turn were used to get a more accurate combined spec- 
trum. This procedure was iterated until the result did 
not change anymore. The resulting spectrum shown in 
Fig. [T] has a S/N ratio of a few tens in the K band and 
a resolution of ~ 4000. The H band spectrum is noisier 
due to shorter observation time and higher extinction. 

The main K band lines (He i 2.112 /im and Br7) are 
clearly identified. But we also detect several weak He i 
fines at 2.149, 2.161 and 2.184 ^m. Similarly, one He i 
and three H i lines are identified for the first time in the 
H band '^. We do not see any trace of He ii 2.189 fj,m. 
With the current S/N ratio, this line would be detected 
if it had a depth of about 1% of the continuum level. 

For the quantitative analysis of this spectrum, we use d 
the atmosphere code CMFGEN ()Hillier k Miller] [l998[) . 
It allows the computation of non-LTE atmosphere mod- 
els including winds and line blanketing and is appropriate 
fo r hot massive stars. A description of th e code is given 
inlHilher fc Miller! (|1998f ). and we refer to iMartins et all 
()2007t l for a summary of the procedure used to build the 
models. To derive the parameters of interest, we have 
run models with: 19000 < T^s < 30000 K, 3.0 < logg < 
4.75, 0.1 < He/H < 2.0 and 10"^ < M < lO'^-^ Mq 
yr~^. These models and synthetic spectra included, in 
addition to H and He, C, N, O, Si, S and Fe. We adopted 



^ These lines were seen bv lEisenhauer et al.l l|2005l ) in the average 
spectrum of several stars, but not in individual spectra. 
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TABLE 1 

Physical properties derived from our quantitative 

analysis for different toff appropriate for, b0-2.5v stars. 
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Fig. 1.— Spectrum of S2 after addition of all the SINFONI 
observations since 2004. The main lines are indicated. 

the solar abundances of iGrevesse fc Sauval ()1998[ ). 

3. THE NATURE OF S2: A GENUINE B STAR 

3.1. Spectroscopic classification 

The spectrum shown in F i g. [T] i s typical of an early 
B star (se e also iGhez et al.| I2003D. Ins pection of the 
atlase s of [Wallace fc Hinklg (|1997D and iHanson et al.1 
(|2005l ) shows that for late O stars, He il 2.189 fim is 
expected. We do not detect this hue in spite of the good 
enough S/N ratio. For stars with spectral types later 
than B3, the He i lines vanish. Since we detect several 
of them, we can safely argue that S2 has a spectral type 
between BO and B2.5. Fig. 12 of IHanson eFaLl ()2005D 
shows that in supergiants and giants, Br7 is clearly sep- 
arated from He I 2.161 iim, while for dwarfs Br7 is broad 
enough to merge with the He i line. This latter morphol- 
ogy is similar to what we observe for S2 (see below for 
a quantitative comparison). Hence, in addition to be- 
ing spectroscopically identified as an early B star, we 
can conclude from these qualitative arguments that S2 
is also a dwarf and not a supergiant/ giant. 

3.2. A mass estimate 

The key question we want to answer here is the exact 
nature of S2. Spectroscopically, S2 is unambiguously an 
early B dwarf. However, this does not necessarilly mean 
that S2 is a young, massive star. It could be the core 
of an older, evolved star which would have lost its en- 
velope through tidal interaction with SgrA* but would 
still spectroscopically look like a B star. In that case, 
the star could have formed far away from SgrA* before 
being dragged to its proximity and experiencing an en- 
velope stri pping. Th ere would be no paradox of youth. 

Recently, |Davies_& King (2005) and lDrav et aLl (|2006f ) 
modeled this stripping process as well as the subsequent 
evolution of the remaining core in the vicinity of SgrA*. 
They showed that, under certain condition on the IMF 
and the capture rate, the population of S stars could be 
explained by tidal stripping. The argument is however 
based on the number of stars predicted at the position of 
the S stars in the HR diagram. Expressed differently, the 
conclusion comes from the fact that the modeled stripped 
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stars can reach the luminosities and effective tempera- 
tures of B stars. But strictly speaking, this does not 
exclude that the S stars are genuine massive stars, which 
would lie at the same position in the HR diagram. 

For that, the only parameter which can unambiguously 
be used is the stellar mass. The most massive stars able 
to experience envelope stripping are AGB stars. Accord- 
ing to stellar evolution, such objects are the evolved de- 
scendents of main sequence stars with M < 8 M© . The 
core of such objects is actually much less massive: Ta- 
ble 4 of iForestini fc Charbonnell (|1997f) shows that their 
mass is lower tha n 1 M©. In the scenario explored by 
iDrav et all (|2006l ). only stars having lost at least 99% 
of their envelope could explain the S stars. Thus, their 
mass should be at most ~ 1 M©. 

To estimate the mass of S2, we can rely on gravity and 
radius: M = gB? /G where g is the gravity, R the ra- 
dius and G the gravitational constant. Gravity can be 
derived from the shape of Br7, while radius is straight- 
forwardly obtained from the knowledge of effective tem- 
perature and luminosity. Below, we explain how we pro- 
ceeded to constrain these parameters. 

Effective temperature. Tcff is usually derived from the 
ratio of lines from two consecutive ionization states: He i 
/ He II lines for O stars and Si in / Si iv for early B stars. 
Unfortunately, S2 does not show He ii and Si lines in the 
K band. Hence, we can only rely on semi-quantitative 
arguments to estimate Teff. We have shown above that 
S2 is a B OV to B2.5V star. Accordin g to the recen t 
studies of iDufton et all (|2006[ ) and lTrundle et all (|2007f l. 
such stars have 19000 < T^s < 30000 K. We have thus 
run models for Toff = 19, 22, 25, 27 and 30 kK. 

Luminosity. To estimate the star's luminosity, for each 
Tcff we adjusted log -^ in our models to match the ab- 
solute K band magnitude of S2. This magnitude is de- 
fined as MK = mK - Ak - DM with mK = 14 
the observed magnitude of S2 (jPaumard et al.l 120061 ). 
Ak = 2.25 the K band extinction at the position of 
S2 (|Schodel et al.l[2007) and DM = 14.50 the distance 
mod ulus for a distance to the Galactic Center of 7.94 
kpc ( Eisenhauer et al.l 120031 ). The resulting absolute K 
magnitude of S2 is -2.75. The range of values we obtain 
for log -j^ is 4.30-4.80. In practice, the uncertainty of 
±0.2 on Ak and 0.5 kpc on the distance to the Galac- 
tic Center translate into an uncertainty of about 0.12 in 

Radius. The radius of the star is simply derived from 
L = A-KCFBR^T'^ff {o'B being the Boltzmann constant). 
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Fig. 2. — Determination of logg for the case T^s = 22000 K. 
The black solid line is the observed Br7 line. The dashed (dot- 
ted, dot-dashed) line is a model with logg = 4.0 (3.00, 4.50). All 
other parameters are kept constant (especially He/H=0.5). The 
wavelength range used for the logg determination is indicated. 

The minimuni radius aUowed by our study (correspond- 
ing also to the minimum mass of the star - see below) 
is estimated for a luminosity reduced by the uncertainty 
(0.12 dex). The values of Rmin are given in Table [TJ 

Gravity, log a was c onstrained from the shape of Br7. 
iRepolust et al.l ()2005f ) showed that its wings (resp. ab- 
sorption core) get broader (resp. weaker) when gravity 
increases. For large log 5, the blue wing merges with the 
He I 2.161/im line as discussed previously. Fig. [2] shows 
how the line profile changes when log (7 increases from 3.0 
to 4.5 for Toff = 22000 K. The best models correspond 
to logg ~ 4.0. A quantification of the goodness of the 
fit in the Br7 region of the spectrum was performed by 
means of a x^ analysis. To minimize contamination by 
Hel, we used only the range 2.164-2.179/z?ti. Fig.[3]shows 
the result of this analysis and confirms that logg — 4.0 
is preferred. From this curve, we also derive a 3 cr lower 
limit to log 5 (log 5mm)- Table [T] shows the results for 
other Tcff. The values o f logg we obtain are all typical of 
a dwarf star {logg^ 4.0. lTrundle et al. 2007 '1. confirming 
our previous "spectroscopic" findings. 

Rotational velocity. We derived a projected rotational 
velocity of 100±30 km s~^. It best accounts for the shape 
of the He i 2.112 /im doublet: for larger V sini the lines 
merge into a single component, for lower values, they are 
too seperated. This value is common ly found for early 
B stars: according to lAbt et all ()2002l ). V sini = 127±8 
(108±8) km s-i for B0-2V (B3-5V) stars. 

Table [T] summarizes, for each T^s, the values of logg 
and radius. Since we want to test if S2 is the core of an 
AGB star, we also give the minimum mass of S2. All our 
estimates are larger than 14.1 Mq. This is larger than 
the initial mass of AGB stars, and consequently is well 
above the mass of the core of such a star. From that, one 
can thus safely conclude that S2 is not a stripped AGB 
star, and is really a genuine B star. 

4. MASS LOSS RATE AND HELIUM CONTENT 
4.1. Mass loss rate 



150 



100 




Fig. 3. — X as a function of logg for the Tcff = 22000 K model. 
The best fit model has logg = 4.0, as also seen from Fig. [2] 

Below a certain threshold, Br7 is insensitive to M. 
When M increases above this limit, the Br7 absorption 
profile is progressively filled by emission. By determin- 
ing this threshold, we can set an upper limit on the mass 
loss rate of S2 (see Fig. H]). In practice, the quantity we 
constrain is the wind density p = M / {AivR^Voo) (where 
R is the stellar radius and Voo the terminal velocity) . De- 
riving a mass loss rates implies we know Voo- Since the 
lines we observe in the S2 K band spectrum are insensi- 
tive to 11(50, we have to assume values. Terminal velocities 
of B dwarfs are difficult to estimate since the winds are 
usually weak and the spectra do not show P-Cygni lines 
or emission lines from which it is usually derived. Hence, 
we have adopted v^ — 1000 km s~^ which is appropriate 
for late O dwarfs (e.g. iBouret et al] |20031. If t;oo would 
be lower than 1000 km s""'^, estimated mass loss rates 
would be lower. Our assumptions thus lead to conserva- 
tive upper limits of M < 3 x 10~^ M© yr^^. This is lower 
than M required bv iLoebl ()2004[ ) to feed SgrA* without 
transport of angular momentum. 

4.2. Helium content 

The near infrared spectrum of 82 is dominated by H 
and He lines. We thus could constrain the He/H abun- 
dance ratio. Qualitatively, when this ratio increases, all 
He I lines have stronger absorption profiles, while H i 
lines get weaker. In practice, we used the He i 2.149 ^m 
and He i 2.184 /ztti lines to constrain He/H since these 
lines are almost insensitive to microturbulence, which is 
not the case of He i 2.112 fjim and He i 2.161 fxm. Fig. 
[5] shows the changes in the He I line profiles when He/H 
is varied between 0.1 and 0.5 (model with Tcfi = 22000 
K). We find that the He/H ratio varies from about 0.45 
at Toff = 25000 K up to as high as 1.20 at T^s = 19000 
K. We proceeded as for logg [y^ analysis) to derive a 3 
a lower limit on the value of He/H. All these lower limits 
are larger than 0.25, confirming that S2 is He enriched. 

With He/H > 0.25, S2 falls into the category of the 
so-called "He rich" s tars, a class of B stars with He/H 
between 0.3 and 10 ()Smith|[l996[ ). Interestingly, these 
peculiar stars have a very narrow distribution of spectral 
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Fig. 4. — Effect of mass loss on Br7. When M increases, Br7 
starts to depart from a pure absorption to an emission profile. 
Characterizing the value of M for which this occurs provides an 
upper limit on the mass loss rate of S2. 

types centered around B2, similar to what we find for S2 
(B0-2.5V). The best studied He-rich star is a OriE. The 
origin of its surface abundance pattern is explained by 
a combination of specific wind properties and magnetic 
field. The star should have a wind weak enough for ion 
decoupling to occur, i.e. the radiative acceleration, es- 
sentially due to metals, is not redistributed among the 
passive plasma (mainly H and He) because of a too low 
densi ty and reduced Coulomb forces (Krticka & Kub|3 
I2OOII) . In that case, Helium might accumulate at the sur- 
face of the star. But this extra Helium can remain on the 
surface only if turbulence, and its associated mixing ef- 
fect, is suppressed. This calls for a magnetic field strong 
enou gh to freeze the stellar surface (Hunger & Grootg 
Il999f) . Our upper limit on the mass loss rate of S2 
is consistent with a relatively weak wind. In this sce- 
nario, chemical inhomogeneities (spots) are expected on 
the surface as a result of the interplay between the mag- 
netic field geometry and the wind. We cannot test this 
prediction with our observations since we do not spatially 
resolve the stellar surface of S2. Similarly, no CNO ab- 
normalities are expected in this model, which we cannot 



test either due to the absence of CNO lines (observed 
and expected) in the near-IR spectrum of S2. 

If the wind is strong enough to develop into a normal 
homogenous outflow, one might still observe a significant 
He enrichment if the s tar has a magnetic fi eld and rotates 
fast enough. Indeed iMaeder fc Mevnetl (200 5) predict 
that a strong (10^ G) magnetic field can lead to sohd 
body rotation of the star. This favors the diffusion of 
species by meridional circulation. As a consequence, a 
star can experience strong sur face He enrichment: Fig. 
10 of IMaeder fc Mevnetl (|2005D shows that a 15 Mq star 
can reach a He mass fraction of 0.3 (around 0.12 in He/H) 
in 12 Myr. Contrary to the previous scenario, large N 
and reduced C abundances are also expected. 

In view of the available models and of the currently 
known stellar properties, we thus tentatively propose 
that S2 is a magnetic star. 
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Fig. 5. — Determination of He/H. The black solid line is the 
observed spectrum of S2. The red dashed (blue dotted) line is a 
model with He/H=0.10 (0.50). T^g = 22000 K and logg= 4.0. The 
He rich model is required to fit the S2 spectrum. 
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